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Abstract 

We study the static potential between electric charges in the finite temperature 
three dimensional compact gauge theory on the lattice. We show that in the de- 
confinement phase at small separations between the charges the potential contains 
a linearly rising piece which goes over into the Coulomb potential as the distance 
between the charges is increased. The linear potential is due to the gas of magnetic 
dipoles which are realized as monopole-anti-monopole bound states. 

1 Introduction 

The compact abelian gauge theory in three space-time dimensions at finite temperature 
has two phases separated by the phase transition [|1], 0, |^. At low temperature the 
electric charges are confined, while at high temperature the confinement disappears. These 
properties of the model originate from the monopole dynamics, the monopoles being 
topological defects appearing due to the compactness of the gauge group. 

In the low temperature phase the abelian monopoles form a plasma of magnetic 
charges. Due to the long-range nature of the gauge fields associated with the monopoles 
the behavior of the plasma is Coulombic. In the weak coupling regime the plasma is 
sufficiently dilute to warrant the use of mean field methods in order to analyse the non- 
perturbative behaviour of the system. As it was shown in Ref. Q test particles of oppo- 
site electric charges being immersed in the magnetic Coulomb plasma experience confining 



forces at large separations. The confinement appears because of the formation of a string- 
hke object between the charges. The string has a finite thickness of the order of the inverse 
Debye mass and a finite energy per unit of string length ( "string tension" ) . 

At high temperature the physics of monopoles changes: monopoles and anti-mono- 
poles form magnetically neutral bound states These states are filling up the vacuum 
with a neutral dipole plasma and the confinement of electric charges disappears. This 
phenomenon can be explained in two ways. First, the field of the magnetic dipoles is 
much weaker at large distances compared to the field of the monopoles and thus the 
dipole is unable to induce a non-zero string tension. The second explanation is that in 
the dilute dipole plasma the Debye screening is absent and effectively this corresponds 
to an infinitely thick confining string. Since the electric flux of the string is a conserved 
quantity, it is no more concentrated in a small region ( "core of the string" ) and thus the 
string "dissolves" as we go from low to high temperatures^. 

Despite the fact that the dipoles themselves are unable to create the confining string, 
they may still be responsible for the non-perturbative physics of electric charges. The 
string tension corresponding to the linear part of the long distance potential of electric 
charges in a pure magnetic monopole gas increases with addition of the magnetic dipoles 
to the monopole gas [§. At short distances the potential of electric charges in a pure 
magnetic dipole gas contains a piece which rises linearly with the distance between the 
charges 

Note that nontrivialities in a short distance potential appear also in the zero-tempe- 



rature theory. According to Ref. |]I0|], at distances R much smaller than the correlation 
length of the monopole-anti-monopole plasma the potential contains a perturbative con- 
tribution due to one-photon exchange plus a non-perturbative piece proportional to i?", 
a ~ 0.6. A non-trivial short-distance potential may have many physically interesting con- 



sequences, see, e.g. Ref. |]1T| for a discussion in the context of QCD and other theories. 
The monopole binding in 3D compact QED is qualitatively similar to the formation of the 
instanton molecules in the high temperature phase of QCD suggested to be responsible 



for the chiral phase transition |T^]. 

The structure of this paper is as follows. In the next section we review some results 
concerning the physics of the dipole gas in the continuum following Ref. Sect. | 
contains the results of our simulation of the compact abelian gauge model in three dimen- 
sions. We show that indeed the short- distance potential contains a linear piece, which 
can be explained as due to the non-trivial dynamics of the gas of lattice dipoles. Our 
conclusions are summarized in the last section. 



^Another mechanism of the phase transition in the Georgi-Glashow model due to magnetic vortex 
dynamics is discussed in Ref. 0. 
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2 Dipole gas in continuum 



A magnetic dipole is a magnetically neutral localized pair of a monopole and an anti- 
monopole separated by a distance r. The magnetic moment of the dipole is Jl = 
where Qm is the magnetic charge of the constituent monopole and r is the relative position 
vector pointing from the anti-monopole to the monopole. If the typical distance between 
dipoles is much larger than the dipole size r, then the dipoles may be treated as point-like 
particles. This condition can be written as follows: 

^ = plr<l, (1) 

where p is the mean density of the dipole gas. We shall see below that in our calculations 
this condition is always fulfilled. 

The action of two interacting point-like dipoles with magnetic moments /i^ and /if, 
located at positions Xa and x^, respectively, is given by the formula: 

Jib] Xa, Xb) = {Pa " 0) {fib ■ 0)0(^^0) (^a " ^b) , (2) 

where 1^(3/)) (x) = l/(47r|x|) is the propagator for a scalar massless particle in three 
dimensions. 

Below we consider the case of the fixed absolute values of the dipole moments, /Iq = 
fi = const. The case of fluctuating dipole moments is considered in Ref. 0. 
The statistical sum of the dilute dipole gas can be written as follows: 



oo /-AT 

^=Y.J^^j d^^i j dVi ■■■ j ^^^N J d^PN exp 



■ ^ N 

-7y H V{fIa,lIb;Xa,Xb) 

a, 6=1 



(3) 



where ( is the dipole fugacity. As in the case of the monopole gas the dipole fugacity 
is a non-perturbative quantity, since ( ~ e~^°, where 5*0 ~ g~'^ is the action of a single 
dipole and ge = 2ti / Qm is the fundamental electric charge in the theory. 
The partition function (^) can be rewritten as follows 0, 

Z = / Dx exp{- / d'x \'-(Bxf - 4.C 2^!^ ] } . (4) 



where \dx\ = \l dx ■ Rescaling the field, x ~^ (^ttC)^^^^ /i^^ X; ^ind the coordinates, 
X — s> {AttQ^^^^x, we immediately realize that the dynamics of the gas is controlled by the 
dimensionless constant 

A = 47rC/x^ (5) 
The vacuum expectation value of the dipole density, Pd(x) = J2iS^^\x — Xj), is 

pa=<p^{x) >=47rC(l + 0(A)), (6) 
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where the last equahty is vahd for small couplings A. In this regime the coupling A is 
proportional to the density of the dipoles, A = Pd/^^, up to O(A^) terms. Therefore the 
smallness of the parameter A can be interpreted as a requirement for the density of the 
dipole moments to be small: 

Pd/i^<l. (7) 

The static interaction of particles with electric charges Qe is the sum of the perturbative 
contribution from the one-photon exchange and the contribution from the dipole gas, 
respectively: 

V{R) = I /^(2D)(i?) + E^^iR) , E^^iR) = - hm 1 < W{Cn.T) > , (8) 

Z 1 — »oo J 

where Crxt stands for the rectangular RxT trajectory of the test particle and D{2D){.R) = 
— (27r)~^ log(mi?) is the two-dimensional propagator for a scalar massless particle, m is a 
regulator of the dimension of mass. We consider here test particles with an electric charge 
equal to one unit of the fundamental charge. The case of particles of arbitrary charges is 
discussed in Ref. 0. 

At small distances between the test electric charges, -C r, the behaviour of the 
potential was shown to be as follows [^]: 

2 

El^{R) = aR[l + 0(R/r)\, a = vr^ ( r = ^p^r , < r . (9) 

where a is given to leading order in A. 

At large distances, R^ r, the classical energy of the electric charges in the dipole gas 
grows logarithmically: 

El^{R) = \og{^) [l + 0{r/R)] , i? » r . (10) 

According to eq. (H) the dipole gas non-perturbatively renormahzes the coupling constant 
Qe at large distances and the full potential (^) has the following behaviour: 

r(i?) = |^D(2D)(i?), e=l + ^A + 0(A2), i?»r, (11) 

where e is the dielectric constant (permittivity). 

The linear term at small distances (^) is a result of the interaction of the dipole clouds 
which surround the external electric sources. Indeed, in three dimensions the static test 
charge trajectories jc may be considered as electric currents running along the "wires" jc- 
These currents induce a magnetic field which encircles the trajectories and is defined to 
leading order by the classical Maxwell equations. The closer to the current the larger the 
magnetic field is. Since the energy of a magnetic dipole becomes lower with increasing 
magnetic field strength, the density of the magnetic dipoles should increase towards the 
position of the electric test charges. Therefore the dipoles form clouds near test particles 
and the interaction of these clouds is responsible for the non-perturbative part of the 
inter-particle potential. 
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3 Compact QED in 3D: Numerical Results 



We consider the 3D compact U{1) gauge model with the standard Wilson action, Sp = 
—f3cos6p, where 6p is the abelian field strength tension constructed from the abelian 
compact fields 6i and f3 is the lattice coupling which is related to the electric charge Qe in 
the following way, 



P 



1 



(12) 



a is the lattice spacing. Our numerical results are obtained on the lattice of size 32^ x 8. 
The finite temperature phase transition corresponds to the lattice coupling Pc = 2.3, 
Ref. 0. The high-temperature phase corresponds to /3 > Pc- 

We do realize that our results can not be expected to work in the confinement phase. 
Therefore, when fits are made only the numerical results for (3 > I3c are used. Still, we 
consider it worthwhile to show the numerical results in a larger domain of (3 values as 
they illustrate the qualitative difference between the two regimes. 

We show the behaviour of the density p of the abelian monopoles vs. (5 in Fig. |l|(a). The 
density p decreases rapidly with increasing of P {i.e., with increasing of the temperature). 
Note that in the deconfinement phase, j3 > (5c the density of the monopoles is non- 
zero. Analysis of the gauge field configurations gives that the monopoles appear as tight 
monopole-anti-monopole bound states in which the constituents are separated basically 
by one or two lattice spacings a. 
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Figure 1: (a) The density of the abehan monopoles, p, vs. lattice coupling [3] (b) vacuum 
dielectric constant (permittivity) e. 



We study the potential between oppositely electrically charged particles measuring the 
correlations of two Polyakov lines, L{x), separated by a distance i?, 

V{R) = - log < L{0)L+{R) > . (13) 



5 



The lattice potential V{R) is fitted by the following formula: 

V{R) = —^D{R) - log|cosh[a {R - 1/2)] | + C , (14) 

where the effective coupling, /3fit, the short string tension, a, and the additive energy 
renormalization, C, are the fitting parameters. The function D{R) = D{R,0) is the 
two-dimensional lattice propagator for a scalar massless particle: 

L—l ^2iTik-x/L 

^^""^ ^ 21? ,^5^3 2 - cos(27rA;i/L) - cos(27rfci/L) ' ^^^^ 

where the zero mode, /ci = /c2 = 0, is excluded. We are using the massless propagator 
since in the high temperature phase (which we are interested in) the potential between 
electric charges is long ranged according to eq. (|ll]) (the dipoles are unable to induce 
a finite correlation length). Function (|1^) fits the numerically obtained potential (|T3|) 



with = x^/d.o.f ~ 1 up to distances ra ~ 8, while at large distances the function 
become unacceptably large. Thus we are fitting the potential (|13|) by the function (|I3|) 
at the distances R = 0, . . . ,8. 

The second term in Eq. (jl^ corresponds to the linear term modified by finite volume 
corrections. In the limit L — > cxd this term transforms to the ordinary linear term aR. 
Note that the finite volume corrections in the first term are already taken into account 
since according to Eq. ([15|), D{R + L) = D{R). 

The coupling /?fit used as a fitting parameter in Eq. ([l^ should be equal to the "bare" 
coupling, P, in first order perturbation theory. However, due to the non-perturbative 
corrections coming from the dipole gas 0, this coupling gets renormalized. The ratio of 
the renormalized and bare couplings, according to Eqs. (|TT|,|T2p is the dielectric constant 
(permittivity): 

e{P) = = , (16) 

P 9e 

which is plotted in Fig. [^(b). In the deconfinement region, j3 > Pc, close to the phase 
transition the dielectric constant e is approximately equal to 1.2. A similar investigation 
of the charge renormalisation in four dimensional zero-temperature QED has been done 
in Refs. |]13| . 



The non-perturbative part of the inter-particle potential is given by the full potential 
minus the perturbative one-photon exchange, 

Vnp{R, (3) = Kxp(i?, P) - -w^DiR) . (17) 

Pfit IP J 

In Eq. (1) the non-perturbative part of the potential is solely due to the dipole gas 
contribution, E^'^lR), and, according to our fits it is linear. We show the full potential 
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Figure 2: (a) Nonperturbative part of the test particle potential, eq. (|T^, vs. full 
potential, for /? = 2.35, and (b) the nonperturbative parts of the potential for various 
values of /3. 



and the non-perturbative part of the potential /5 = 2.35 in Fig. 0(a). Indeed, we observe 
that the non-perturbative part is linear and is substantially smaller than the perturbative 
one. The non-perturbative parts of the potential for various (5 are shown in Fig. @(b). 

The tension of the short string a, defined in eq. (|l^ is plotted vs. P in Fig. |^(a). 
Note that a is a decreasing function of P (or, equivalently, of temperature). The short 
distance string tension at the high temperature side is non-zero due to the magnetic 
dipoles dynamics discussed in Sect. ^, Eq. (^. 

In order to check the consistency of the results obtained, with the dipole gas picture 
we first estimate the dipole size from the string tension and the monopole density using 
eq. (|^). Taking into account the fact that the density of the dipoles, pd, is half the density 
of the constituent monopoles, p, we get: 

8CT 



We evaluate r by this formula and plot it as a function of j3 in Fig. ^(b). In the decon- 
finement phaseQ, /3 > Pc, the dipole sizes become smaller as /3 increases, in qualitative 
agreement with the theoretical estimates. Note that at large P the monopole sizes are of 
the order of 4 or 5 lattice spacings while the observed sizes are of the order of 1, . . . , 2 
spacings. 

There are two reasons which may explain observed quantitative difference in r of a 
factor of 3. First, the distribution of the dipole magnetic moment for a real dipole state is 
unknown while in our theoretical estimations we have assumed fixed dipole moments for 
definiteness. Another reason might be that the analytical formulae for the dipole density 



^In the confinement phase formula ( [l8|) should not work since in this phase a fraction of the monopoles 
is not bounded in dipole states. 
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Figure 3: (a) the average distance between constituent monopoles in a dipole state, 
eq. (|18D; (b) the short distance string tension obtained with the help of the fit (|l^). 



d^) and the string tension (P) are not apphcable. This might be caused by a violation 
of the dipole diluteness condition (|I|) or of the requirement for the density of the dipole 
moments to be small, eq. (^. 

The quantity characterizing the diluteness of the gas, ^, eq. (|ip, is of the order of 
0.1 for the measured distances between monopole constituents. However, the dipole sizes 
found with the help of eq. ([18|) give ^ ~ 0.3. The last relation indicates that the point- 
like dipole gas approach might be "on the edge of applicability". However, the density 
of the dipoles as well as the dipole size decrease with the temperature and therefore the 
condition (|I]) should be satisfied better as temperatures increases. 

The other assumption which has been used in the derivation of the analytic expressions 
for the short string tension a, eq. (|^), and the permittivity e, eq. (pH]), is the smallness 
of the parameter A, eq. (^. Using Eqs. we can express A as a function of the 

monopole density, the string tension and the lattice coupling: 

A^5i^. (19) 

With the help of this equation we get that A is a decreasing function of P which is equal 
to 0.8 3X (3 = Pc and decreases to 0.3 at /5 = 2.5. Thus near the phase transition and 
even at the largest studied values of (3 the requirement fails to be fulfilled. However, 
at large (3 the coupling A is about to be small enough for the theoretical formulae for the 
short string tension a, eq. (^, and the permittivity e, eq. (p!TD, to start working. 

Thus we may conclude that our theoretical formulae may work qualitatively rather 
than quantitatively. As we have seen above this is indeed the case: the linear potential is 
observed while the mismatch in the self-consistency check for the dipole sizes for the largest 
value of /? is a factor of 3. Better agreement between theoretical and numerical values 



can be observed for the dielectric permittivity. Indeed, according to eq. (|TTD , — 1 = A/ 



3 ~ 0.118 while numerical data shown in Fig. [I|(b) gives the result e„„m — 1 = 0.177. 
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4 Conclusions and Acknowledgments 



We have observed the nonperturbative piece in the inter-particle potential in high tem- 
perature compact QED in three dimensions to rise linearly with the distance. This effect 
is in qualitative agreement with the predictions of the pointlike magnetic dipole gas model 
with fixed dipole sizes. The electric permittivity of the vacuum was also calculated and 
turns our to be very close to the theoretical value. 

These results may have interesting applications for the physics of the gauge theo- 
ries possessing monopole topological excitations which form bound states. One of these 
theories is the electroweak model in which the formation of the Nambu monopole-anti- 



monopole pairs has been observed in the high temperature phase |14 
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was partially supported by the grants RFBR 99-01-01230a and CRDF award RPl-2103. 
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